Quantum Phase Transition in Extended Attractive Hubbard Model 
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We have made a variational analysis on an evolution of superconductivity from weak to strong 
coupling regime. In contrast to a crossover without thermodynamic anomaly found in a dilute 
system, we show the existence of a quantum phase transition near half filling. The transition is 
driven by charge density waves instabilities. We have found that superconductivity and charge 
density waves coexist in the presence of a weak intersite repulsion. The ground state phase diagram 
is determined and the quantum phase transition is an attractive version of the Mott transition. 



The Landau Fermi liquid theory Q explains thermody- 
namic properties and responses to external fields of inter- 
acting fermion systems, and predicts the existence of col- 
lective excitations. It is based on the assumption of adi- 
abatic continuity from the noninteracting system. How- 
ever, this assumption is not always satisfied in strongly 
interacting electron systems. On increasing repulsion, a 
metallic Fermi liquid (FL) phase breaks down, and the 
system becomes an antiferromagnetic (AF) insulator [Q. 
This is the Mott transition — a kind of quantum phase 
transition. We investigate whether a corresponding tran- 
sition exists in a superconducting (SC) phase on increas- 
ing attraction instead of repulsion. 

The Bardeen Cooper Schrieffer (BCS) theory has been 
the basis for our understanding of superconductivity and 
is known to be valid in the weak coupling regime 
We think that the BCS theory in the strong coupling 
regime should be critically examined in relation with the 
small coherence length superconductors discovered in the 
last decades In the strong coupHng case, tightly 

bound pairs of electrons are expected to behave like Bose 
particles called bipolarons, and Bose-Einstein condensa- 
tion (BEC) should take place. According to the work 
of Nozieres and Schmitt-Rink |^ the evolution with in- 
creasing attraction is a smooth crossover. This result 
is widely accepted, at least in qualitative sense. The 
smooth evolution scenario is justified in a dilute system, 
but must break down near half filling (n = 1) where 
a charge density waves (CDW) phase is expected. We 
reexamine the evolution taking the possible occurrence 
of CDW phase into account, and find a quantum phase 
transition between SC and CDW. Moreover, we find that 
these different kinds of orders compete, and coexist un- 
der certain conditions. A coexistence of SC and CDW 
has been recently observed in a stripe phase of high Tc 
superconductors 

We employ an infinite dimensional (d — oo) extended 
attractive Hubbard model on an AB-bipartite hypercu- 
bic lattice l^-Q to simplify the problem, without losing 
essential feature of quantum fluctuations. 
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The last term represents a nearest neighbor repulsion 
which plays an important role to stabilizing CDW. 

To improve the BCS theory in a controlled way, we 
rely on a variational method as a guiding principle for 
the ground state. We consider the following Gutzwiller 
mean-field variational wave function (GMF) 
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where |MF) is a mean-field (MF) wave function. This 
is an extension of the Gebhard wave function to in- 
clude SC order with |MF). Parameters ^ia-, rji, and are 
determined so as to eliminate any self-energy diagrams 
in d = oo. The Gutzwiller projection incorporates local 
electronic correlations, and > 1 in an attractive sys- 
tem. We have found significant deviations from the BCS 
theory with increasing U, and g^— >2asC/— >oo |16|. 

Before discussing SC we comment on the normal state. 
We substitute the noninteracting Fermi sea (FS) for 
|MF), and the resulting |GMF) is the Gutzwiller wave 
function |GWF). On increasing U, the discontinuity q at 
the Fermi surface of the momentum distribution function 
diminishes. This means that the quasi-particles effec- 
tive mass increases. At U = Ubk = Seq and arbitrary 
filling n, with eq — \/2/7rt* the average band width, 
a quantum phase transition arises. We obtain q ^ 
for U > Ubr. The electrons completely localize form- 
ing bipolarons, where the spin susceptibility Xs — > and 
a spin gap opens up. This quantum phase transition 
corresponds to the Brinkman-Rice metal-insulator tran- 
sition in the repulsive case [0. It does not actually 
occur if we take SC into account. However, it implies 
that a crossover from FL to a Bose liquid takes place 
at [/ ~ Ubr- This interpretation is consistent with the 
peak positions of superconducting condensation energies, 
as shown in Fig. ^ 

For SC we substitute in the BCS wave function |BCS) 
for |MF) in Eq. (|), and we caU it the Gutzwiller BCS 
wave function |GBCS). We obtain an optimized varia- 
tional energy with the s-wave symmetry m& . We show, 
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in Fig. 0, SC condensation energies /S.E = i^pL — Esc of 
the attractive Hubbard model {V* = 0). AE increases 
following the BCS formula ~ t* exp{—at*/U) in the weak 
coupling regime, where a is of the order unity, decreases 
proportional tot* ^/C/ in the strong coupling regime, and 
has a peak near U ^ Ubr independent of n. Simple 
|BCS) and |FS) cannot reproduce this behavior, and we 
must take Fermi liquid effects into account for both FL 
and SC. Therefore, the work of Nozieres and Schmitt- 
Rink Q is insufficient due to a lack of Fermi liquid effects. 
Note that this by no means suggests a break-down of the 
BCS theory. In the original BCS theory, FL is consid- 
ered as the normal state, and shows that FL is unstable 
towards SC in the presence of the infinitesimally weak 
attraction between quasi-particles Then, operators 
in the BCS-reduced-Hamiltonian, should be considered 
as dressed quasi-particles. Instead, we start from the at- 
tractive Hubbard model represented by bare electrons, 
and incorporate Fermi liquid effects by the Gutzwiller 
projection. We have shown that an infinitesimally weak 
attraction actually leads to SC by |GBCS). AE is a rem- 
iniscent of an expected SC transition temperature Tc in 
a crossover from BCS to BEC |,|jl|jl|. 

We show in Fig. ^, the SC order parameter |(cixCi|)| 
and the energy difference between |BCS) and |GBCS). 
Deviations are most prominent in the intermediate cou- 
pling regime. SC order parameters are overestimated by 
|BCS) compared with |GBCS). In |BCS), variational en- 
ergies are lowered only by the presence of an SC order 
parameter, whereas in |GBCS), we can further lower vari- 
ational energies by Fermi liquid effects. The agreement in 
a weak coupling limit is expected, but the fact that they 
also agree in a strong coupling limit, is a result of d = oo 
limit. In a strong coupling limit U oo, the attrac- 
tive Hubbard model maps onto the spin model 111, i.e. 
the Heisenberg model with an external magnetic field. 
However, in d = oo, nearest neighbor spi n-sp in interac- 
tion reduces to only Hartree contributions , and a MF 
treatment becomes exact. We have checked correct MF 
limits, (nij-riii) — > n/2 and |(ci|Ci-|-)| ^ n{2 — n)/2, as 
[7 — > oo. In finite dimensions d < oo, the simple spin 
waves theory {1/S expansion) is sufficient to give lower 
order parameters compared with a MF theory, and we ex- 
pect that deviations persist finite as J7 ^ oo. Although 
|GBCS) shows a significant deviation from |BCS), as far 
as a pure SC is concerned, the evolution turns out to be 
continuous, which is consistent with common wisdom. 

Now, we explore a discontinuous evolution. SC and 
CDW are characterized by the diagonal and off-diagonal 
long range order in the Nambu density matrix 



|SC + CDW) = Y[[ 
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and we investigate a competition of these orders. In or- 
der to take a coexistence phase (SC+CDW) into account, 
we put 



as |MF), where ako- denotes an annihilation operator for 
CDW bogolons defined by 
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and Q — (tt, tt, • • •) is a commensurate CDW wave vec- 
tor. Our wave function |GMF) comes from a natural uni- 
fication of the spin-bag approach ||2^ and the resonating- 
valence-bond idea |^l| for the attractive system. Optimal 
functional forms of variational parameters u and v are an- 
alytically determined beyond a MF theory, and further 
minimizations are performed numerically ]l^ . 

In the absence of V* , CDW and SC are degenerate at 
n = 1, and minimum variational energies do not change 
for a coexistence phase. This result was obtained in the 
MF theory [Q and is not altered by the Gutzwiller pro- 
jection. This reflects a symmetry |^ of the attractive 
Hubbard model. However, such a macroscopic degener- 
acy is lifted under any weak perturbation. In fact, for 
V* > we obtain a pure CDW ground state at rt = 1, 
and we find a coexistence phase at n ^ 1. 

We show, in Fig. ^ order parameters and variational 
energies as functions of n. CDW suppresses SC but these 
orders can coexist near n — 1. Doped holes in CDW are 
expected to behave like Bose particles in a strong cou- 
pling regime. As far as concentrations of holes are di- 
lute, they can not completely destroy CDW. Moreover, 
a dilute Bose gas exhibits BEC at zero temperature, and 
as a result the coexistence is realized. Further reduction 
of n induces a quantum melting of CDW, and we find a 
pure SC. The quantum phase transition is of the second 
order. 

A coexistence of SC and charge stripe orders was ob- 
served by NQR measurements in cuprates Although 
dimensionalities and CDW structures are quite different 
from the experiments, our results qualitatively agree with 
observed behaviors of order parameters. The coexistence 
of d-wave SC and AF has also been found by the vari- 
ational Monte Carlo method in the two dimensional t-J 
model 1^^. Neutron scattering experiments of cuprates 
revealed that three orders (SC, CDW, and AF) coex- 
ist near 1/8 hole doping Therefore, our results are 
complementary to previous theoretical results and are 
supported by experimental observations. A coexistence 
phase is also observed in bismuthate superconductors by 
infrared measurements [p4[ . However, in bismuthates, SC 
is realized after several tens of percent of doping jlj,^ . 
To understand these experiments, we must take into ac- 
count a finite dimensionality, incommensurate CDW, lo- 
calization effects due to impurities, and so on. 

Another possibility for the ground state is a phase 
separation 123] between electrons rich and poor regions. 
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In fact, the inverse of the charge susceptibility = 
dii/dn shows that the system is a proximity of it. 
However, unlike neutral bosons, charged electrons are 
hard to exhibit a macroscopic phase separation due to 
the long range Coulomb interaction. 

We show in Fig. ^, phase diagrams at fixed values of U 
or n. Clearly, CDW are more stable near n = 1. The MF 
theory underestimates CDW instability. A pure CDW is 
stable at n = 1 and V* > 0. We find a quantum phase 
transition between a pure SC and a coexistence phase 
by varying U and/or V*. |MF) gives a reentrant behav- 
ior of the phase boundary. In contrast, |GMF) yields 
smaller value of V* for CDW with increasing U. Our re- 
sults can be understood by a simple picture: on increas- 
ing U coherence length of Cooper pairs becomes smaller, 
leading to a smaller value of V* for CDW due to larger 
susceptibility of different pairs. This behavior is also ex- 
pected from an effective Hamiltonian in the strong cou- 
pling regime, the anisotropic Heisenberg model Q: SC 
and CDW orders correspond to ferromagnetism in the xy 
plain and AF along the z direction respectively, and the 
effective exchange interactions are Jxy = 2t*^/Ud and 
= 2t*'^ /Ud + V* /d. Therefore, an infinitesimally weak 
V^* > is sufficient to break a pure SC in [/ — > oo. 

The phase diagram of the infinite dimensional ex- 
tended attractive Hubbard model (0) is shown in Fig. 
|. At = ill or C/ = [|6|, our result is the same 
as the MF results. The rest of the phase boundaries are 
different from those of simple MF theories. The fact that 
we have found a quantum phase transition from a pure 
SC, proves the existence of a discontinuous evolution. 
The BCS theory is one of the most successful theories in 
physics, so this is a significant phenomenon. The quan- 
tum phase transition is an attractive version of the Mott 
transition. 

In summary, we have made a variational analysis on an 
extended attractive Hubbard model in d = cx). We have 
found the importance of electronic correlations through 
evaluations of superconducting condensation energies. 
As far as a pure superconducting phase is concerned, 
an evolution with increasing attraction turns out to be 
continuous in a dilute system. However, such a smooth 
evolution scenario breaks down near half filling, where 
charge density waves take place. We have found a quan- 
tum phase transition between a pure superconducting 
phase and a coexistence phase of superconductivity and 
charge density waves induced by varying interactions 
and/or filling. This phenomenon can be viewed as a Mott 
transition in an attractive system. 

We appreciate enlightening discussions with Prof. I. 
Terasaki. This work is supported by Waseda University 
Grant for Special Research Project (98A-855, 99A-564). 
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FIG. 1. Superconducting condensation energies Ai 
= i?pL — -Esc of the attractive Hubbard model at n — 1 
0.75, 0.5, 0.25, compared with several differences at n = ] 
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FIG. 2. Improvements on |BCS). SC order parameter 
calculated with |GBCS) are shown ai n — 1, 0.75, 0.5, 0.2f 
combined with the BCS results (the dotted line). The inse- 
is the energy difference. 
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FIG. 3. Coexistence of SC and CDW at [/ = 3eo and 
V* — O.Ieo. A coexistence phase is denoted by SC+CDW, 
and pure phases are denoted by SC or CDW. The upper Fig. 
shows order parameters of CDW Sriu = (uact ~ nBcr)/^ (the 
solid line) and SC |ci|Ci|| (the dotted line, A and B refer to 
the sublattice). The lower Fig. shows superconducting con- 
densation energies. 
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FIG. 4. Quantum phase transition between weak and 
strong coupling superconductivity. The solid line is the phase 
boundary of a pure SC and a coexistence phase (SC+CDW). 
The dotted line is the MF result. The left Fig. is obtained at 
17 = 5eo, lOeo. At n = 1, a pure CDW is stable for V > 0. 
The right Fig. is obtained at n = 0.75, 0.5. Normal states 
are stable at (7 = 0, and there exists a transition between a 
homogeneous FL and a pure CDW. 




FIG. 5. Phase diagram of the extended attractive Hub- 
bard model in d = CO. The curve is the phase boundary of a 
pure SC and a coexistence phase (SC+CDW). Italic indicates 
a phase on a plain. 
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